Modification of radiation pressure due to cooperative scattering of light 
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Cooperative spontaneous emission of a single photon from a cloud of TV atoms modifies 
substantially the radiation pressure exerted by a far-detuned laser beam exciting the atoms. 
On one hand, the force induced by photon absorption depends on the collective decay rate of 
the excited atomic state. On the other hand, directional spontaneous emission counteracts 
the recoil induced by the absorption. We derive an analytical expression for the radiation 
pressure in steady-state. For a smooth extended atomic distribution we show that the 
radiation pressure depends on the atom number via cooperative scattering and that, for 
certain atom numbers, it can be suppressed or enhanced. Cooperative scattering of light by 
extended atomic clouds can become important in the presence of quasi-resonant light and 
could be addressed in many cold atoms experiments. 

I. INTRODUCTION 

Coherent radiation by large collections of small particles has been studied extensively with 
applications in antennas, atmospheric scattering or plasma diagnostics. The behavior of those 
macroscopic scattering devices is ruled by the interference of the fields scattered by many micro- 
scopic elements. In atomic physics, matter-light interactions have been studied extensively in the 
past decades, ranging from spectroscopy to laser manipulation of the atomic internal and exter- 
nal degrees of freedom. Intriguing properties can be obtained when many atoms are interacting 
with photons, as studied in the seminal work by Dicke [l| and more recently by cold samples of 
bosonic atoms in a trap 0]. Other complex quantum features are obtained when the quantum 
properties of light and/or of atoms are taken into account, as is done e.g. in studies of quantum 
cryptography and quantum computation 01 • Often a quantum description of two level atoms and a 
quantum description of the light field in terms of photon operator is very convenient to describe the 
atom-light interactions. However one should not forget that many features can well be described 
with classical models for the atoms with a polarizability and for the light field. Many features of 
Dicke superradiance can thus be explained using classical antenna theory. As atoms appear to be 
excellent systems to study any possible deviation from classical many body features, it is a general 
interest to understand and to control cooperative effects, even at a classical level, in a cloud of 
atoms. New intriguing effects can arise, when fluctuations due to the coupling with the vacuum 
modes can no longer be described by a classical field approach and the atoms can also become 
entangled during the cooperative scattering. Such cooperative scattering with atomic ensembles 
can appear in a number of experimental situations in free space [1] or in cavities [H]. 

The problem of collective spontaneous emission has recently received growing interest with the 
study of single photon superradiance from N two-level atoms prepared by the absorption of one 
photon of wave vector ko = &o e z [6-9]. It has been shown that the photon is spontaneously emitted 
in the same direction of the incident photon with a cooperative decay rate proportional to N and 
inversely proportional to the size of the atomic cloud [id ]. These studies considered the decay of 
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atoms prepared in the 'timed Dicke state': 
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where \g%, g 2 , ■ ■ ■ , &?> ■ ■ ■ ,9n) is a Fock state in which the atom j is prepared in the excited state 
e and all the other atoms being in the ground state g, and Tj is the position of the atom j. In 
these papers, only resonant photons have been considered and their reabsorption by other atoms 
has been neglected. As a situation closer to an experimental realization, we consider a cold atomic 
cloud with a Gaussian density distribution irradiated by a plane wave laser beam (see Fig. [T]) . If 
the laser beam is detuned far enough from the atomic resonance to avoid multi-atom absorption, 
a timed Dicke state is generated. The large detuning also ensures that the dominant fraction of 
photons is elastically scattered. 

In this paper, going beyond the mere consideration of the light emission by an atomic ensemble, 
we will revisit the phenomenon of radiation pressure in a regime where scattering is dominated 
by cooperative effects. Mechanical effects combined to cooperative emission have been studied for 
atomic ensembles in cavities [H| or in optical lattices However, they also occur in free space 
[lH ] and, as we will show in this paper, can modify not only the emission properties but also the 
absorption process. Indeed, the absorption properties are affected by superradiant enhancement of 
the spontaneous decay rate. On the other hand, cooperative scattering bundles the scattered light 
into forward direction, which leads to a suppression of radiation pressure. This suppression is only 
limited by the inhomogeneity of the atomic cloud, which occupies a finite volume. For volumes 
larger than the Dicke limit, |rj| 3> A, the net effect of cooperativity is a suppression of radiation 
pressure far below the value expected for independent scattering. 



FIG. 1: (Color online) Scheme for collective scattering of a single photon from a pump laser beam ko by a 
spatial distribution of atoms. 

We proceed as follows. We solve the Schrodinger equation for an ensemble of N two- level atoms 
with the resonance frequency u a driven by a far-detuned laser field with the frequency uq = oj a +Ao 
and the wave vector ko. For a laser field sufficiently detuned from resonance so that the average 
number of excited atoms in the cloud is less than unity, t/NQq < Ao (where Qo is the Rabi 
frequency of the laser field), we assume that only one among N atoms is excited by absorption of 
a single laser photon, which is then spontaneously emitted into vacuum modes with frequencies 
LOk = uj a + Ak and wave vectors k. Within the Markov approximation (valid for t ^> R/c, where 
R is the dimension of the cloud [101 ]) we obtain expressions for the time-evolution of the state 
amplitudes allowing us to derive a generalized expression for the steady-state radiation pressure in 
the presence of cooperative effects. 

In the limit where the atomic cloud can be described by a smooth spherical Gaussian den- 
sity distribution, an analytical calculation of the structure factors allows us to derive an explicit 
expression for the radiation pressure as a function of volume, shape, atom numbers, and pump 
laser detuning. In particular, we identify two regimes distinct by different scalings of the radiation 
pressure with atom numbers. 
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II. TIME EVOLUTION OF THE SYSTEM 

The atom-field interaction Hamiltonian is in the rotating-wave approximation (RWA) [3] 

N ro N 
H = hJ2 -fa,e iAot - iko ^ +h.c. +hY,Y, [5fc^4 eiAfc *" ik ' r3 ' + h - c -] • ( 2 ) 
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Here, H.q is the Rabi frequency of the interaction between an atom and the pump mode (which is 
treated as a classical field), &j is the lowering operator for atom j, is the photon annihilation 

operator, and gf, = d^Jio/(heoV p h) describes the coupling between the atom and the vacuum modes 
with volume V p h- The j-th atom has lower and upper states denoted by \gj) and | e j ) , respectively. 
With this approach we treat the atoms as simple two-level systems. We also assume that all 
atoms are driven by the unperturbed incident laser beam, thus neglecting dephasing by atoms 
along the laser path or by near field effects, which could arise for large spatial densities. Calling 
|0) a = \9i,--,9n) the atomic ground state and \j) a = |<?i, .., ej, .., gjv) the state where only the atom 
j is excited, we assume that the total state of the system has the following form: 

N 

|tt(t)> =a(i)|0) a |0) k + £ /3 J (t)|j)a|0) k + ^ 7 kW|0)a|l)k • (3) 

j=i k 

The above expression assumes that only states with at most one atomic excitation contribute to 
the effects described in this paper. 

The time evolution of the amplitudes is obtained by inserting the Hamiltonian ([2]) and the 
ansatz © into the Schrodinger equation, dt\^(t)) = —(i/H)H\^(t)): 

N 

a(t) = - i -n e iA ° t Y,me~ ik °' Ij , (4) 

1 3=1 

fat) = -*-n oa (t)e- iA ° t+i ^ - *£^7k(*)e- 4Z M-H k - r , > (5) 
1 k 

N 

7kW = -^e iAfci E^( i ) e " ik ' rj • ( 6 ) 
3=1 

Integrating Eq. © over time and substituting 7k (i) in Eq. ([5]) we obtain N coupled equations: 

N c t 

fat) = --n a(t)e- lAot+i ^ -Y,9kY. e ik - (rj " rm) / e^^^CO^ • ( 7 ) 
1 k m=l Jo 

We have numerically checked that due to the presence of the driving term, the solution quickly 
evolves toward a "driven timed Dicke state" 0], characterized by 

BJt) = ^e- iAot+iko ^ . (8) 



Once inserted the ansatz ([8]) in Eqs. dl]),([7]) and going to continuous momentum space via J2k 
V ph (2-7r)- 3 J dk, we get 



(9) 
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$(t) = - l -^NQ oa (t) + iA p(t) - ^ J dk g 2 k \S N (k)\ 2 x jT* e^K-**)*^* _ ^ j (10) 

where <SV(k) = jj Ylj=i e _ *( k ~ k °)" r -? is the structure factor of the atomic cloud. We can now apply 
the Markov approximation and write 
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Then, Eq. (jlOp becomes 



$(t) = - l -VNQ a(t) + (*A - ~TN SN \ p(t) , (12) 

using the definition of T = (V^/Trc)^^ . Here, using spherical coordinates, we introduced the 
quantity sjy: 

s N = ^- / 27F # I* d6 sin6 \S N (k , 6, (f))\ 2 . (13) 

47T Jo Jo 

In steady state and neglecting saturation, yiVOo ^ A , we may assume a(t) ~ 1 and find the 

P - « ^ . (14) 

III. FORCES IN THE MARKOV APPROXIMATION 

The two terms in the Hamiltonian ([2]) yield two different contributions to radiation pressure 
force: 

F aj + F ej = -V rj H . (15) 

We will be interested in the average absorption force, F a = J2j(^aj), and emission force, F e = 
jj J2j(Fej), acting on the center of mass of the whole cloud, F a + F e = m&cM, where slcm is the 
center-of-mass acceleration and m the mass of one atom. 

The first term, F a j = ^Tiko^lol&je 1 ^ 01-1 ^ rj — h.c] results from the recoil received upon absorp- 
tion of a photon from the pump laser and has an expectation value on the timed Dicke state given 
by: 

F a = (F aj ) = ^^Im[a(t)P*(t)) . (16) 



The second contribution, F e j = iJ2k.^9k[^jd\ L e lAkt lkr J — h.c], results from the emission of a 
photon into any direction k. The expectation value on the general state ([3]) is: 

(F ej )=^a%[ft(t) 7 J( f )e iA ^-c.c.] . (17) 

k 

Substituting the time integral of 7k(£) from Eq. ([6]) and inserting the timed Dicke state from Eq. (|8|) 
we obtain the average emission force: 

F e = -E k hkg 2 k \S N (k)\ 2 x \fi(t) ; * e^-°)''/3*(t - t')dt' + c.c.j (18) 
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In the Markov approximation and going to continuous momentum space, 

F e =-\m\ 2 ^Jdk(hk)gl\S N (k))\ 2 S(k-k ) (19) 
Denning /n = j^ Jq 2 " d(j) J$ dOsinOcosO |Sjv(&0i <t ) )\ 2 1 we find 

F e = -hk r\(3(t)\ 2 f N . (20) 

Finally, using Eq. and a « 1 in Eqs. (fl"6]) and (f20j) . the steady-state value of the total radiation 
pressure force on the center of mass of the atomic cloud is 

F c = F a + F e = hkoT 4A *™° T 's% - f N ) ■ (21) 

Both forces, F a and F e , are directed along the direction of the incident laser beam. 

This is the main result of the paper. We identify contributions from cooperative absorption 

and cooperative emission. The common prefactor can be obtained from the standard single-atom 

n 2 

radiation pressure force F± = hkoT ^$^_ r2 , by substituting in the fraction the natural linewidth 
by the collective linewidth, T — > NTsn, and the Rabi frequency by the collective Rabi frequency, 
fio ~~ * v^V^o- Additionally, the cooperative radiation pressure force is weighted by the difference 
of structure factors, sn — fx, where the sn part corresponds to the cooperative absorption process 
and the /at part to the cooperative emission. We will see that both forces may cancel completely 
in some cases. 



IV. STRUCTURE FACTORS FOR SMOOTH DENSITY DISTRIBUTIONS 

Assuming a smooth Gaussian density distribution with ellipsoidal shape, n(r) = noexp[— (x 2 + 
y 2 )/2a 2 — z 2 /2a 2 ], we can evaluate the structure factor by replacing the sum J2j by an integral 
/ dr n(r), and obtain 
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(22) 



where a = ka r and rj = a z ja r is the aspect ratio. Setting k ~ ko, one can show that for elongated 
clouds, rj > 1, 



Aa^f] 2 - 1 
fin) 

J OO 



erf 



a{2if - 1) 



erf 



1 



a 



r/ 2 - 1 S °° 4a 2 1 6 



(23) 



For spherical clouds (rj = 1) one finds 



/o ° " 4<r 2 



Aa 2 
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(24) 



For cr, r\ 3> 1, can be approximated by Soo ~ s 00 yfirFe t * [1 — erf(F)], where F = a/rj = ka 2 /a z 

(w) 

is the Fresnel number. For large Fresnel numbers Soo — > Soq. 
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FIG. 2: (Color online) The red solid line shows the iV- dependence of the radiation pressure force resulting 
from pump photon absorption only, for r\ = 1, a = 5 and A = 10 3 T. The green solid line traces the 
radiation pressure force, when photon emission is taken into account, and the blue dashed line traces the 
single atom radiation pressure. 



V. DISCUSSION 

For a smooth-density spherical cloud the ratio between the cooperative radiation pressure force 
(|2ip and the single-atom force becomes 

F ° - 4A o + r2 Ns U-h>\ (25) 
Fx ~ 4A 2 + iV 2 r 2 s^ V sj • (25) 

The emission pattern is described by the term. An isotropic emission pattern, obtained e.g. for 
small volumes, leads to vanishing f^. In contrast, for large volumes the recoil at emission com- 
pensates the recoil at absorption, ~ Sqo, which results in mainly forward emission. The above 
expression also shows that even for isotropic emission pattern, cooperative effects can strongly 
modify the radiation pressure by a modified cooperative absorption (see Fig. [2]). The dependency 
of Eq. (|25h on atom number exhibits a maximum at N cr = 2Ao/rsoo = 4cx 2 (2Ao/r). Assuming 
large detuning An S> T, the modification of the radiation pressure force for small atom numbers, 
N <C N cr , is 

F c 2N 

On the other hand, for large atom numbers, N 3> N cr , the radiation pressure becomes independent 
of the cloud size, 

The maximum radiation pressure, obtained at N = N cr , is F c jF\ ~ An/2r<7 2 and hence, for large 
volumes with a > \J An/2r, below the radiation pressure expected for uncorrelated scattering. For 
small volumes the cooperative can exceed the uncorrelated radiation pressure. 

It is interesting to compare this with the Dicke limit of small clouds, <C A, where the 
structure factor becomes 5jv(k) ~ 1. Consequently, the surface- integrated values are sjy = 1, 
/at = and N sr = 2An/T, so that the force ratio for small and large atom numbers becomes, 
respectively, 

Yr N and T^n^—) • (28) 

Note that the phase shift of the laser beam after passing through the atomic cloud can affect 
the driving term of Eq. ([2]). Indeed, for a spherical atomic cloud with resonant optical density 
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60 = (3A 2 /2-7r) J dz n(0, 0, z) = 3N/a 2 , the characteristic atom number N cr can be expressed in 
terms of the phase shift experienced by the pump laser beam on its path across the cloud: A(j> = 
6or/(4Ao) = 6N/N cr . In other words, for atom numbers so large, that despite its large detuning 
the pump laser beam is shifted by more than 2ir, the radiation pressure force is dramatically 
suppressed. For low spatial densities of the atomic cloud, this phase shift remains small compared 
to the free propagation phase shift ka z and only weakly modifies the scattering properties of the 
atomic cloud. As we assume that each atom is only driven by the external field, we also neglect 
multiple scattering of a photon in the atomic cloud. The above calculations are thus not valid for 

b > 4Ag/r 2 . 

The collective radiation pressure described in this paper is likely to be of importance for samples 
with large on-resonance optical thickness (bo 3> 1). Samples with such an optical thickness are e.g. 
used to study multiple scattering of light in cold atoms [15] (with bo « 40) and most Bose-Einstein 
condensates also fulfill this criterium. For an experimental observation of the collective radiation 
pressure, a low intensity quasi-resonant laser would be required to measure radiation pressure, 
e.g. in a time of flight experiment after the exchange of several photons per atom on average. 
Such experiments are state of the art as well for cold atom in magneto-optical traps as for atoms 
from a Bose-Einstein condensate. For Bose-Einstein condensates the modification of the scattering 
properties of light might also allow to address the problem of spurious heating in optical dipole 
traps where spontaneous scattering might appear in some limits. 

We note again that cooperativity in the radiation pressure as described in this pap er does not 
rely on self-organization or correlations between different scattering events (as in ll|, |12|]) but it 
is intrinsic in the scattering of a single photon by a correlated system of N atoms. Although the 
calculated expression of the average radiation pressure does not exhibit specific quantum features 
and could be in principle obtainable classically, the study of higher moments of the atomic distri- 
bution and fluctuations should shed light on the quantum properties of the entangled state [8, 16]. 
This topic will be at the aim of our forthcoming studies. 



VI. CONCLUSION 



We have shown by an explicit calculation that the radiation pressure depends on the atom 
number via cooperative scattering. The assumption of a smooth density distribution allowed us to 
integrate the structure factors. We have shown that the radiation pressure is strongly affected by 
both a cooperative modification of the absorption and the emission properties of the atomic cloud. 
We have identified different regimes, where enhanced absorption can increase the radiation pressure 
for intermediate atom numbers or where forward emission almost cancels the recoil imparted to the 
atoms at the absorption. These calculations indicate that experiments are in reach with state of 
the art clouds of cold atoms. As the driven timed Dicke states involved in our model are entangled 
atomic states, radiation pressure might become an interesting tool to investigate non classical 
features of cooperative scattering of light by clouds of cold atoms. 
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